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Abstract. We show that the governing equations for two-dimensional grav- 
, ity water waves with constant non-zero vorticity have a nearly-Hamiltonian 

i structure, which becomes Hamiltonian for steady waves. 
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J> . 1. Introduction 

The mathematical study of water waves was initiated within the framework of 
linear theory with the work of Airy, Stokes, and their contemporaries in the nine- 
teenth century. Periodic two-dimensional water waves are of special interest since 
the typical water waves propagating on the surface of the sea (or on a river or lake) 
present these features. Stokes [23] noticed that actual water wave characteristics 
deviate significantly from the predictions of linear theory. This started an extensive 
i-£h ■ study of the nonlinear governing equations for water waves. 

A celebrated development in water-wave theory was the discovery by Zakharov 
' [29] that the governing equations for two-dimensional irrotational gravity water 

waves have a Hamiltonian structure - see the discussions in [2, 13, 14]. The aim 
°f this paper is to present a nearly-Hamiltonian formulation for two-dimensional 
gravity water waves with constant vorticity. For irrotational flows (zero vorticity) 
we recover Zakharov's result. Moreover, if we restrict our attention to steady 
waves, a Hamiltonian structure emerges; we refer to [17] for an in-depth discussion 
5_i ■ of the Hamiltonian structure of steady irrotational water waves. While Zakharov's 

ideas were generalized in various directions by several authors - see [5, 19] for a 
survey of contributions in this direction, the elegance and simplicity of the nearly- 
Hamiltonian formulation for flows of constant vorticity makes it conceivable that it 
might be instrumental in deriving qualitative results for such flows. Related to this 
possibility, recently [12] stability results for steady water waves with vorticities that 
depend montonically on the depth were derived from the variational formulation 
provided in [11]. The reason for considering water waves with constant vorticity is 
twofold. Firstly, by the Kelvin circulation theorem [18, 20] a two-dimensional water 
flow that is initially of constant vorticity remains so at later times and therefore the 
restriction to such flows is justifiable. Secondly, as already pointed out, the elegance 
of the nearly-Hamiltonian formulation valid within this setting is mathematically 
attractive. From the physical point of view we notice that while irrotational flows 
are appropriate for waves propagating into a region of still water [20] , water waves 
with vorticity describe wave-current interactions - see the discussions in [6, 24, 25]. 
Tidal flows are the most prominent example of water flows with constant vorticity 
[15]. 
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2. Preliminaries 

To describe two-dimensional periodic water waves it suffices to consider a cross 
section of the flow that is perpendicular to the crest line. Choose Cartesian coor- 
dinates (x, y) with the y-axis pointing vertically upwards and the x-axis being the 
direction of wave propagation. Let (u(t, x, y), v(t, x, y)) be the velocity field of the 
flow, let {y = 0} be the flat bed, and let {y = r](t,x)} be the water's free surface. 

For gravity water waves the restoring force acting on the water's free surface is 
gravity and the effects of surface tension are neglected. Assuming the water density 
to be constant (p = 1) - this is physically reasonable cf. [20], we obtain the equation 
of mass conservation 

(1) U X +Vy=0. 

Appropriate for gravity waves is also the assumption of inviscid flow [20], so that 
the equation of motion is Euler's equation 

(2) 



U t + UU X + VUy = -P x , 

v t + uv x + Wy = -Py - g, 



where P(t, x, y) is the pressure and g is the gravitational constant of acceleration. 
The free surface decouples the motion of the water from that of the air so that 

(3) P = Patm on y = t](t,x), 

must hold, where P a tm is the atmospheric pressure [18]. Since the same particles 
always form the free surface, we have 

(4) v = T) t + wq x on y = r](t, x). 
On the flat bed the boundary condition 

(5) v = on y = 0, 

expresses the fact that water cannot permeate the rigid bed y = 0. The governing 
equations for periodic two-dimensional gravity water waves propagating over a flat 
bed are (l)-(5), with the specification that the periodicity is reflected in the fact 
that all functions u, v, P, rj exhibit a periodic dependence in the x-variable of, say, 
period L > 0. Other than the nonlinear character of the equations, the main 
difficulty in their analysis lies in the fact that we deal with a free-boundary value 
problem: the free surface y = rj{t,x) is not known a priori. Throughout this paper 
we consider flows of constant vorticity, that is, throughout the fluid the vorticity 

(6) lo = v x - Uy 
is constant (u> e R) throughout the fluid domain 

n(t) = {(x, y) e M 2 : < x < L, < y < r)(t, x)}, 
the free surface of which is given by the graph 

S(t) = {(x, y) e K 2 : < x < L, y = rj(t, x)} 
restricted to a period cell. Furthermore, we require that 

(7) / u(t,x,0)dx = 0, t>0. 

Jo 

The relevance of this last condition is explained in Section 3. 
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3. The nearly-Hamiltonian formulation 

For two-dimensional flows the incompressibility condition (1) ensures the exis- 
tence of a stream function ip(t, x, y) determined up to an additive term that depends 
solely on time by 

(8) U = tPy, V = ~1P X . 

Since (5) becomes ip x (t,x,0) = 0, we can determine tp uniquely by setting ip = 
on the flat bed y = 0, that is, we set 

rVa 

(9) ip(t,x Q ,y Q ) = u(t,x ,y)dy for (x , y ) G il(t). 



This explicit formula shows that ip is x-periodic with period L. In terms of the 
stream function, the vorticity u> is determined from (6) by 

(10) Aip = -uj in Q{t). 

Let us now introduce the (generalized) velocity potential <p(t, x, y) via 

(11) u = tp x -ujy, v = ip y . 

Notice that this is not the standard reduction by the Weyl-Hodgc decomposition 
(see [4, 28]) since the vector field W = (—uiy, 0) is divergence free but does not 
satisfy the boundary condition W ■ n = on a free surface that is not flat (n being 
the outward normal to the boundary). The value of the (generalized) potential tp at 
(t,Xo,yo) can be determined by integrating ip x along the horizontal segment with 
endpoints (0, 0) and (xo, 0), and subsequently tp y along the vertical segment joining 
(x ,0) to (x ,y ) G Q(t): 

rxa nuo 

<p(t,x ,yo) = <p(t,0,0)+ (p x (t,x,0)dx + / <p y {t,x ,y) dy 

Jo Jo 

rx r-ya 

= ^(i,0,0)+ / u(t,x,0) dx + / v(t,x Q ,y) dy. 
Jo Jo 

The potential <p is a harmonic function since 

(12) Aip = in Sl(t). 

Notice that, independently of the additive time-dependent term up to which ip is 
uniquely determined by (11), we have by the above explicit formula for <p that 

pxo+L 

ip(t,x + L,y ) - <p(t,x ,yo) = u(t,x,0)dx, (%o,yo) G Cl(t). 

J x 

In view of the x-periodicity of the function u, we see that the right-hand side equals 

L 

u(t, x, 0) dx. At this point some results on steady water waves (water waves for 

o 

which the free surface r], the pressure P and the velocity field (u, v) exhibit an 
(x, t)-dependence of the form x — ct, where c ^ is the wave speed - that is, in 
a frame moving at speed c these waves are stationary) are of relevance. In the 
irrotational case, the existence of steady waves of large amplitude (that is, waves 

that are not small perturbations of a flat surface) satisfying / u(t, x, 0) dx = 

Jo 

was established in [1] - see also the discussions in [7, 26]. For these waves, called 
Stokes waves, we therefore have that the velocity potential ip is i-periodic in the 
cc-variable. On the other hand, for flows of constant vorticity w G R, there are 
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steady waves of large amplitude for which / u(t, x, 0) dx > cf. [9, 10]. For 



these waves therefore the (generalized) velocity potential ip is not periodic in the in- 
variable. Thus, while ip x and (p y are both ^-periodic with period L, the (generalized) 
potential ip is not necessarily L-periodic in the x-variable. The above discussion 
shows that (7) is the necessary and sufficient condition for p to be L-periodic in 
the x-variable. For steady waves, the relation (7) means that the wave speed is 
defined as the mean velocity in the moving frame of reference in which the wave 
is stationary. Imposing (7) for the large-amplitude steady waves with vorticity 
studied in [10, 8, 11, 12, 16, 27] simply means that along the continuum of waves 
constructed in [10] the wave speed is not fixed a priori, but varies according to 
(7). As stated in the Introduction, throughout this paper we consider only flows 
satisfying (7) so that <p is L-periodic in the x-variable. 

In terms of the functions <p and ip, we can recast the Euler equation (2) in the 
form 

V ^t + JlVVf + P + uip + gy 



0. 



Thus at each instant t > 0, 



(13) 



ft + ^iwf + p + + 



gy 



is constant throughout fl(t). 



This is the generalization for flows of constant vorticity of Bernoulli's law [22] for 
irrotational flows (co = 0). In view of (3), we deduce that 



p t + ^\Vip\ 2 +cuip + gr] 



is constant on the free surface S(t). 



Since p is uniquely determined by (11) up to an arbitrary additive term that is 
solely time-dependent, we use this freedom to absorb into the definition of ip a 
suitable time-dependent term so that 



(14) <p t + - |VV>r+tW + <W = on S(t). 

Let £ denote the evaluation of ip at the free surface, 



(15) 



£(t,x) = ip(t,x,r)(t,x)), t>0, x e [0,L]. 



Given the constant vorticity oj, the functions £ and rj, taken to be smooth and L- 
periodic in the X-variable, completely determine the motion. Indeed, the function 
x i ^ r){t,x) fixes the fluid domain Cl(t), and £(f, •) is the appropriate boundary 
data for a linear elliptic problem of mixed type that determines ip(t, •, •) at time 
t > 0. At any given time t > 0, fix rj and £, and let ip be the unique solution of the 
boundary- value problem 



(16) 



Aip = in Q(i), 
<p = £ on S(t), 
y y = on y = 0. 

that is L-periodic in the x-variable. Knowing ip, we determine the velocity field 
(u, v) from (11), the stream function ip from (10), and the pressure P from Bcrnoidli's 
law (13). 
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The total energy of the wave motion in a period cell is given by 

(17) H=[[ {^±^+ gy }dydx. 

In the above expression, the first term represents the kinetic energy (energy of 
motion), while the second term is the gravitational potential energy (energy of 
position). We now show that H is completely determined by the functions £ and 
r]. In order to do this, we introduce the Hilbert transform, T(rj), defined as follows. 
At any given time t > 0, fix r\ and £, and let ip be the unique solution of the 
boundary-value problem (16). Since (8) and (11) ensure that {ip + ^J/ 2 ) is the 
harmonic conjugate to ip, the Hilbert transform T(i]) £ of £ is given by 

(18) (nv)t){x) =x(t,x) + ^ V 2 (t,x), xe[0,L], 
where 

(19) X (t,x)=i>(t,x,ri(t,x)), t>0,xe[0,L]. 

The equations (18)- (19) show that x is completely determined by rj and £ via 

(20) x = T(r 1 )Z-^r i 2 . 
Let 



1 

n 



be the outward normal to the free surface y = T](t, x). With dl = yjl + r\ x dx, using 
Green's identity, we get 



H = \ff \V^\ 2 dydx+ l - f grfdx 



/fi(t) 

(21) =7,1 V' [-tpxVx + i> y ] dx + ^ If ipdydx+\ f grfdx 

z Js(t) z J Jsi(t) z Jo 

in view of (10) and the fact that ip — on y = 0. Another application of Green's 
identity for the functions ip and y 2 yields 



J J ^ ( " w y 2 - 2V>) dydx = J ^ (y 2 Aip - ip A(y 2 )) dydx 

= J & In ~ ~in^) dl = j (-Vxtpx + ipy) V 2 dx - 2 xvdx 



>s(t) 

if we take into account (10) and (19). Therefore 

(22) // ipdydx = ~l rfdx-\l (-n x ip x + ip y ) i] 2 dx + / xvdx. 
J JQ(t) "Jo 1 Jo Jo 

Taking into account (8), (11), we further transform the expression (21) for H into 

^J s ^(4>-7;y 2 )[ i PyVx + Px-uy}dx + ^J^ (g - ^- rj) rj 2 dx + ^ xvdx 
in view of (22). Since 

6x = Vx + <PyVx, X G [0, L], 
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and using (18), we obtain 

rL -, i-L 



H(,,,() = 1/ t, •T(i,)£<fe+ i | 9>, 2 <ti 

-V' 

6 7o 



r-i 2 /•£ 



(23) -| jf i, I/" f/.r + ~~ / // ! 

We now present the main result of this paper. 

Theorem 1 The governing equations for periodic two-dimensional gravity water 
waves of constant vorticity lu are equivalent to the nearly- Hamiltonian system 

5H 



(24) 



SH 



with y = rj being the free surface, £ being the evaluation of the (generalized) velocity 
potential on the free surface, H = H(rj,£) being the total energy of the motion, and 



being the evaluation of the stream function on the free surface. 

Remark For irrotational flows (u> — 0) the system (24) is Hamiltonian: we 
recover Zakharov's result [29]. □ 

Before proceeding with the proof of Theorem 1, let us review briefly the con- 
cept of an infinite-dimensional Hamiltonian system - see also [14, 21]. An infinite- 
dimensional Hamiltonian system is a system of partial differential equations of the 
form 

(25) f t = J grad H(f), 

where f(t) describes a path in a Hilbert space endowed with an inner product (•, •}, 
the Hamiltonian function H : V — > R being defined on a dense subspace of the 
Hilbert space, and with J being a skew-adjoint operator. The gradient in (25) is 
taken with respect to the inner product (•, •} on the Hilbert space. If the operator 
J is invertible, this set-up defines a symplectic structure on the Hilbert space. The 
system (24) is nearly-Hamiltonian: one can view u as a parameter measuring the 

deviation of (24) from a Hamiltonian structure of the form (25) with / = (^j on 

the Hilbert space L 2 [Q, L] x L 2 [0, L], with operator 

1 N 



1 ' -1 0, 

the Hamiltonian function H being given by (23) with 77 and £ in the dense subspace 
V C L 2 [0, L] of smooth L-periodic functions. 

SH 

Proof of Theorem 1. Let us first compute the ^-gradient of H, — — . We vary £ 
and keep r\ fixed. If 6 is a harmonic function in ft(t), L-periodic in the x- variable 
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and with y — on y = 0, let 60 be the evaluation of 9 on S(t). If is the harmonic 
conjugate of 9 with ^ = on y = 0, from (18) we get 

T(V) [£] = X + \ V 2 , T(V) lZ + e9 }= X + e(+^r, 2 , 
denoting by £ the evaluation of \& on Therefore 



( ^ o) = 1^0 i 



(26) 



~\ J {** x ^ + ^ Xx ~ W 71 Vx ^ 6 ° ^ dx 



in view of (23) and using the periodicity. Since 

6 X = ty y , 9 y = throughout tt(t), 

we obtain 



/ ( x £dx = (V x + V y r] x )<pdx= (-9 y + 9 x i] x )tpdx 
Jo Js(t) Js(t) 

(27) =-/ — ipdl = - 9 — dl= 9 (if x r] x - <p y ) dx 

Js(t) dn J s(t) On J s(t) 



SH 1 f 

(-ff>0o) = —z {<PxVx - <f y +Xx - wrm x )9 dx = / {v-u^ x )9 , 

<>? 1 J Sit) Jo 



/S(t) un JS(t) un JS(t) 

using in the next to last step Green's identity for the harmonic functions <p and 9, 
which satisfy 9 y = <p y = on y = 0. From (26)-(27) we infer that 

r-L 

, , .1 dx. 

>S(t 

if we take into account (8), (11) and (19). Thus 

SH 

— = V - UT} X . 

Since r\ = r] t , we can recast (4) as 

. SH 

SH 

To recover the remaining part of the system (24) , we have to compute -j— . To 

o-q 

do this, instead of working with the expression (23), it is simpler to compute the 
variation in (17). This calculation is not straigthforward due to the fact that there 
is an implicit nonlinear dependence of ip upon 77. As an example of the intricacies 
of the calculation, notice that, contrary to a possible first impression, the Hilbcrt 
transform T(-q) is not skew-adjoint if the surface y — i](t, x) is not flat (see [13, 14] 
for a discussion). 
Let 



(28) 



£i(t,x) = ip x (t,x,ri(t,x)), 
60,2/) = <Py(t,x,ri(t,x)), 



be the evaluations of the partial derivatives of the potential ip on the free surface 
y = r](t, x). Notice that £1 is not to be confused with the function £ x since 
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In terms of the potential <p, we can rewrite (17) as 

L r\ L r) 



H(V,¥>) = ~ J J VVxdydx + ^ J J uj 2 y 2 dydx 



L i) L rj 

(29) + \J J\Vp\ 2 dy+±J J grfdx. 



We now vary the function r\ describing the free surface by 6r), keeping £ fixed. Since 
p and (tp + — y) are harmonic conjugate functions, by analytic continuation [3] the 

function ip has a harmonic extension across the boundary. Varying the domain Q(t) 
to fl e (t) by keeping •) fixed means that we solve instead of (16) the problem 

'Atp e = in Q £ (t), 
p e =i on S e (t), 
dyifi £ = on y = 0, 

where 

Ct £ (t) = {(x,y) e K 2 : < x < L, < y < rj(t,x) + s (6r))(t,x)}, 
S e (t) = {(i,!/)eR 2 :0<x<i,pt|((,x)l£({#,4 

Since 

(p e (t,x,t] + eSr]) = ip(t,x,rf) = £{t,x) 

we deduce that 

ip E (t,x,rj) - <p(t, x,rj) p e (t,x,r]) -<p e (t 7 x,r] + e5r]) 

e e 

= - Sr] ■ d y ip £ (t,x,Tf) + 0(e) 

— > — 6r] ■ <p y (t, x, rj) as e — > 0. 

Therefore, if [Sip]* is the evaluation of the variation Sip of ip on y — r](t,x), we 
showed that 

(30) [<M # = - & fy- 
Other than a variation (5<p of (/?, we have the variation 

l v L L 

SH = — J J lo y 5(p x dydx — J wr/d Srjdx + i ^ uj 2 r] 2 5i]dx 



L t) L L 

(31) + y J (Vp) ■ VStp dydx + 1 y (£ 2 + e|) <fy cfe + J grjSrjdx 



o o 

of H, since 5V<£ = V(5<p. Using the formula 



^ y F(x,y)dy = J F x (x,y) dy + F[x,rj\ rj x 
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the first term on the right-hand side of (31) can be rewritten as 



l r, 



J j ' uySip x 



dydx = 



o o 



toy Sip dy 



dx 



J \jr]{5tp}* r) x 



dx 



(32) 



J wrj[6tp]#r} x 



dx = u> Tfflx&Sv dx 



by (30) and periodicity. Combining (31)-(32), we obtain 

L 



dx 



L v 



(33) 



(Vip) • VSip dydx . 



o o 



Since ip is harmonic in fl(t), applying Green's identity to the pair of functions ip 
and Sip, we infer that the last integral in the above expression equals 

/ Sl P IT" dl = ~ I 6 [6 - £i Vx] Sv dx 
in view of (5) and (11), respectively (28) and (30). Therefore (33) becomes 

rL . 1 e 2 _l C 2 . 



SH 
Thus 



(34) 



£ 2 + £ : 



2 _i_ tf2 



^77 ~ '""^ ~ ' 2 ' ' 2 
Combining (14) with (8), (11), (19), and (28), we obtain 

£ 2 + £ 2 1 
so that (34) becomes 

(35) — = -ip t - ur)Vx& - il + r]x€i& - 

Differentiating (15) with respect to t, we get 

& = + &Vt on y = r](t, x), 
in view of (28). Furthermore, (4), (11) and (28) yield 

£2 = % + (£1 - wr?)ri x on y = r/(i, a;). 
From the previous two relations and (35) we infer 

SH 

*T = " 6 " ^ 

which completes the proof. 



□ 
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4. HAMILTONIAN STRUCTURE FOR STEADY WAVES 



In this section we consider the case of steady water waves with constant vorticity. 
If c 7^ is the speed of the wave, (4) becomes 

-ip x = -cr] x + ip y r] x on y = n{x - ct), 



if we use (8). Thus d x 



so that in this setting 



ip(x — ct, r/(x — ct)j — cr)(x — ct) 

the function [ijj — cy] is constant on the free surface. On the other hand, using (1), 
wc deduce that 



j\u-c) 
Jo 



dy 



= (u(x - ct,r](x - ct)) - cj rj x {x - ct) + J u x dy 
= (u(x — ct, n(x — ct)) — cj rj x (x — ct) — J v y dy 
= (u(x — ct, r)(x — ct)) — cj 7] x {x — ct) — v(x — ct, r)(x — ct)) 



I" 

Jo 



has to equal zero by (1). Since ip y = u throughout the fluid and ip = on the flat 
bed, we have 

rv 

(u -c)dy=(ip- cy) 

10 y=r)(x-ct) 

The constant value k of the expression on the above left-hand side is the relative 
mass flux of the flow [10] . Since field evidence and laboratory measurements indicate 
that for steady water waves that are not near the breaking state the relation u < c 
holds throughout the fluid [20], we have that k < 0. These considerations show 
in view of Theorem 1 that the governing equations for steady water waves with 
constant vorticity lo and wave speed c are equivalent to the nearly-Hamiltonian 
system 

5H 



(36) 



SH_ 

Si] 



— ui(k + ct]) 



where k < is the relative mass flux of the flow. From (36) we readily obtain the 
Hamiltonian formulation for steady water waves with constant vorticity. 

Theorem 2 The governing equations for steady L-periodic two-dimensional 
gravity water waves with constant vorticity u are equivalent to the Hamiltonian 
system 

5_H_ 



(37) 



V 



. = _5H 
Si] 



Here c ^ is the speed of the wave, y = rj is the free surface above the flat bed 
y = 0, ti, is the evaluation of the (generalized) velocity potential on the free surface, 
and 



H = H + uok 



n dx 



coj r o 

Wo " 



dx, 



where H — H{rj,£) is the total energy of the motion. 
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Remark For the Hamiltonian system (37) we have T{rj)^ — en — —rf 2 = k u\ 
view of (20), since £ = k + en. □ 
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